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Abstract 

We analyze the relativistic dynamical properties of Keplerian and non-Keplerian circular 
orbits in a general axisymmetric and stationary gravitational field, and discuss the implications 
for the stability of co- and counter-rotating accretion disks and tori surrounding a spinning 
black hole. Close to the horizon there are orbital peculiarities which can seem counterintu- 
itive, but are elucidated by formulating the dynamics in terms of the orbital velocity actually 
measured by a local, zero-angular-momentum observer. 
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1 Introduction 



Accretion onto supermassive black holes has emerged as the paradigmatic model for the energy 
source of AGNs and quasars (Lynden-Bell 1978, Rees 1984), a picture recently given compelling 
support by HST spectrography of the nucleus of M87 (Harms et al. 1994). 

Detailed development and understanding of such models requires insight into the subtle fluid 
and plasma dynamics close to the edge of a spinning black hole. For near-maximal spin, accretion 
disks can extend almost down to the horizon (modulo the possibility of still poorly understood bar- 
like instabilities, see Blaes 1987). In this strong-gravity regime, general-relativistic corrections can 
produce non-Newtonian behaviour which, at first sight, is strange and counterintuitive (see Allen 
1990). For instance, angular momentum is transported inwards rather than outwards by a viscous 
disk (Anderson and Lemos 1978); there is a reversal of the Rayleigh criterion for stability (which 
conventionally requires the specific angular momentum to increase outwards) (Abramowicz and 
Prasanna 1990); and spinning balls of fluid undergoing slow collapse actually become more spherical 
(rather than more flattened) in the latter stages of contraction toward black holes (Chandrasekhar 
and Miller 1974). 

An insightful series of papers by Abramowicz and co-workers (Abramowicz and Lasota 1986, 
Abramowicz 1990, Abramowicz and Prasanna 1990) reveal how these seemingly disparate anoma- 
lies can be understood in terms of a common, non-Newtonian mechanism: an effective reversal 
of centrifugal force inside a region called the rotosphere, which lies within radius r = 3m for a 
Schwarzschild black hole. (This line of thought is further pursued and expounded in Abramowicz 
1992, 1993 and Abramowicz and Szukiewicz 1993.) 

Viewing a thing from more than one angle often adds depth and flexibility to one's under- 
standing. The group of anomalies that fall under the heading "Abramowicz effect" can also be 
"explained" without overturning familiar Newtonian preconceptions (de Felice 1991, 1994; Page 
1993). Indeed, one very simple explanation, which we shall present in this article, actually exploits 
the close formal resemblance between the Newtonian and general-relativistic descriptions (when 
suitably formulated). The basic remark is that the general relativistic formulae become simplest 
and almost indistinguishable from Newtonian when expressed in terms of quantities that a local 
observer actually measures. 

The gist of the matter can be explained in a few words. The outward rocket thrust needed to 
hold a spacecraft in a circular orbit of radius r around a central mass is proportional to {y\ — v 2 )/r 
according to Newtonian theory, where v and vk are the actual and Keplerian orbital velocities. This 
remains true in Einstein's theory, for a spherisymmetric field, if v, vk are interpreted as velocities 
measured by a local stationary observer. 

This at once makes it clear why, within r = 3m (the orbital radius of a circling photon in the 
Schwarzschild field), the thrust needed is always outward, no matter how fast the spacecraft orbits, 
because here the Keplerian velocity is tachyonic: \vk\ > c. This is the essence of the Abramowicz 
effect. 

Moreover, this way of looking at the phenomenon carries over trivially to the axisymmetric 
case when the central mass is in (steady) rotation. (By contrast, on other interpretations this 
extension is less straightforward (Abramowicz, Nurowski and Wex 1993).) The thrust needed to 
stay on an equatorial circle is now proportional to (v + — v)(\v_\ + v), where v + and t>_ are the 
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prograde and retrograde Keplerian velocities. All velocities are understood to be measured by 
a local observer ("ZAMO") orbiting with zero angular momentum. Because of rotational frame- 
dragging, the two Keplerian velocities are not equal and opposite; they become lightlike at different 
radii, the retrograde orbit at the larger radius. 

Such effects of frame- dragging on counter-rotating orbits have an interest which may be more 
than merely one of principle. Recent observations and numerical simulations suggest that counter- 
rotation may be a less unusual feature in old galactic nuclei than was thought just a few years 
ago (Binney and Tremaine 1987). The gas stream produced by a retrograde encounter between 
galaxies can evolve over billions of years into a counter-rotating disk of stars superimposed upon 
and streaming collisionlessly through the old disk. NGC4550 and NGC7217 are examples of such 
"two-way galaxies" (Rubin et al. 1992, Kuijken 1993). 

In Sec. 2 we introduce the basic concepts and techniques in the simple context of a spherically 
symmetric gravitational field. This introductory discussion clears the way for a concise presentation 
of the corresponding results for equatorial circular orbits in a general stationary axisymmetric field 
in Sec. 3. When the mass of the accretion disk or torus is small compared to the central black hole, 
the gravitational field is well approximated by the Kerr metric. In this case, explicit results are 
available, and are summarized in Sec. 4. Further developments and details are assembled in two 
Appendices. 

2 Spherisymmetric fields 

We begin by considering orbits in a general static, spheri symmetric geometry, described by the 
line-element 

ds 2 = V 2 {r) dr 2 + r 2 dtt 2 - V 2 {r) dt 2 . 

Without essential loss of generality we may assume the orbit confined to the equatorial plane 
6 = 7r/2. We write u = if = dip/dr for the angular velocity as measured by a comoring observer, 
using proper time r. 

Keplerian orbits are timelike geodesies, i.e., extremals of the action / Ldr, where 

L = V 2 i 2 - r V - rfr 2 

and L = 1 on the extremal curve. For a circular Keplerian orbit of radius r, the angular velocity 
uk is thus found from dL/dr = to be given by 

rvk V'(r) m 
1 + rW K V ' 1 ' 

With respect to a local stationary observer, a unit test mass on a circular orbit has 3-momentum 
(i.e., spatial projection of 4-momentum j9 M = dx^/dr) 

P = (9< PV )*P V = ru 

and 3-velocity 

v = p/(l+p 2 )^ = ru/(l + r 2 w 2 )* (2) 
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in relativistic units (c = 1). 
Hence (1) yields 

vl/r = V'(r)/V(r) (3) 

for the Keplerian 3-velocity vk- (In the Newtonian limit V ~ 1 + V^ cwt /c 2 , (1) and (3) reduce to 
the expected familiar results.) 

The orbital radius r = r p h of a freely circulating photon (|i>.k-| = 1) is determined by 

l/r ph = V'(r ph )/V(r ph ). 

In the case of the Schwarzschild geometry V 2 = 1 — 2m/ r and r ph = 3m. Inside this radius, 
Keplerian velocities are tachyonic (v 2 K > 1), corresponding to spacelike geodesic trajectories, dr^ < 
0. 

We now turn to general (non-Keplerian) circular orbits. The mechanical force or rocket thrust 
needed to hold a unit test mass stationary in such an orbit is given by F M = Sp^/Sr (the absolute 
derivative of 4- momentum), a purely spatial and purely radial vector. Its radial component is 

(d L (dL\_dL\_l d_L 
~5t*~ 2 [9rr) \dr\dr) dr J ~ 2 V dr' 



This yields for the (signed) magnitude of (outward) force 

F — - ~ rjvV ~ ' \V 
Recalling (1), this can be re-expressed as 

uj 2 



= (9rr^F r = ^-{l-co 2 r(^--r)\. (4) 




(5) 



K , 



Setting uj = 0, we obtain the force F stat = V'/rjV needed to hold up a stationary unit mass. 
At the photon radius r = r p h, we have uj 2 k = (d(p/dr) 2 K = oo, hence F = F stat is finite for every 
timelike circular orbit and velocity-independent (Abramowicz and Lasota 1986). 

The subtracted term in (5), proportional to u 2 , might be thought of as "centrifugal force." 
Were we to adopt this terminology, we would at once make contact with Abramowicz 1990. In the 
rotosphere, Keplerian circular orbits are spacelike, hence oo 2 K < in (5) and "centrifugal force is 
reversed." 

But there is an alternative way to write (5). Using (2) to eliminate uj in (4) in favour of v, and 
recalling (3), gives 

F=-^4- (6) 
r/r 1 — v A 

In this form, "centrifugal force" (now most naturally interpreted as the subtracted term v 2 /i]r(l — 
v 2 ) in (6)) never changes direction, since (1 — v 2 ) is positive for every physical (timelike) orbit. 
Nevertheless, (6) shows as before that to orbit in the rotosphere (where v\ > 1) always requires an 
outward thrust F, no matter how fast one moves. Indeed, speeding up is self-defeating: the thrust 
needed actually increases with speed in the rotosphere and becomes infinite at the speed of light. 
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Within the rotosphere, it takes less outward thrust to hold a body stationary than to maintain it 
in orbit. One might say that the special- relativistic increase of "mass" with velocity (given by the 
denominator of (6)) literally outweighs the effect of centrifugal force in the rotosphere (Page 1993, 
de Felice 1991). 

It is this property that underlies the reversal of the Rayleigh criterion for stability in the ro- 
tosphere. Consider a toroidal distribution of material (e.g. incompressible fluid) held in steady, 
differential rotation around a gravitating mass by a prescribed non-gravitational force field. Sup- 
pose the specific angular momentum of the distribution increases outwards from the axis. An 
orbiting ring of material, if displaced inwards axisymmetrically, will preserve its angular momen- 
tum and will therefore orbit faster than its new surroundings. In the rotosphere, the outward force 
that would be called for to support it in this new orbit is therefore necessarily larger than that which 
actually supports the surrounding material, and this is as much as the local force field provides. 
The displaced ring is forced to sink further, and this triggers an instability. 

The alternative "explanations" of Abramowicz and Page, however helpful and suggestive, are, 
of course, at bottom merely different forms of words to clothe the formulae. One may adopt either 
(or neither), depending on circumstances and personal taste. One advantage of Page's form, as we 
shall now see, is that it is easy to adapt (6) (but less easily (5)) to the case where the central mass 
is rotating. 



3 Stationary axisymmetric fields 

We pass now to orbits in a stationary axisymmetric asymptotically flat geometry. Let 

£ ( Q t) = dx a /dt, ^ = dx a /dip 

be the timelike and axial Killing vectors (assumed to commute). An observer with 4- velocity 
u a = dx a /dr has angular momentum 

t = p ip = u- f (<p) = g^u . 

For a free-falling (geodesic) observer, this is conserved. 

Any observer orbiting in an azimuthal circle has a 4-velocity of the form 

«* = ir'fcg, + n*&))i (7) 

where Q = dip/ dt is his angular velocity as measured by a stationary observer at infinity. The 
orbiting observer is a zero angular momentum (ZAM) observer (Bardeen 1970 a, 1973, Thorne, 
Price & MacDonald 1986) if u • = 0, i.e., if he orbits with the Bardeen angular velocity 

£i B = -(f(v>)-f(t))/(£(vo-£(v)) = -g<pt/g<p<p- 

The normalizing factor in (7) is then 

UzAM = V=(-g tt )-v. 
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The line-element is now conveniently expressible as 

ds 2 = g AB dx A dx B + p 2 sin 2 6 (dtp - Vl B dt) 2 - V 2 dt 2 (8) 

where x A = (x 1 ^ 2 ) = (r,0) and all metric coefficients depend on r and 9 only. 

We assume that this geometry also has equatorial symmetry, and we shall be interested in 
circular orbits in the equatorial plane 6 = ir/2. As in the previous section, their properties are 
derivable from the Lagrangian 

L = V 2 i 2 -p 2 u 2 - rfr 2 
in which the dot denotes differentiation with respect to proper time r, 

uj = ip — Qst 

is the locally measured angular velocity relative to the local ZAM observer and all coefficients are 
functions of r only. The 3- velocity v of a circularly orbiting particle as measured by the local ZAM 
observer is 

v = pcu/(l + p 2 uj 2 )z. 

Keplerian circular orbits satisfy dL/dr — 0, L — 1, r = 0. This yields a quadratic equation with 
roots 

v+ = -y/tf + P), v_ = -i/{5-(3) (9) 
giving the 3-velocities of prograde and retrograde Keplerian orbits. We have defined 

a = d r \np, f3 = -\(p/V)d r Vt B ., 7 = <9 r lnV, 

5 =( ai + (3 2 ) 1 2. (10) 

The retrograde Keplerian orbit becomes lightlike when 

v- = -1 7 - a = -2/3, 5=ha + -f). (11) 



For the prograde orbit, 



v + = +l 7 -a = 2/3, 5 = ^(a + 7 ). (12) 



The thrust required to hold unit mass in a non-Keplerian circular orbit is given by the absolute 
derivative = 5 2 x^/5r 2 . The only nonvanishing component is calculable from 

F r = 5 2 r/5r 2 = \rf 2 dL/dr. 

This gives for the (signed) magnitude of outward thrust F = (g rr )^F r : 

F = _a(v-v + )(v-v.)_ 

T] 1 — V 2 
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Equations (9) and (13) are invariant, as they should be, under reparametrization of the arbitrary 
radial co-ordinate r. A more obviously geometrical form emerges if we introduce the element of 
proper radial distance 

ds r = (g rr ) 5 dr. (14) 

Then 

F = j > (15) 

p 1 — v 2 ds r 

The conical deficit factor dp/ds r (slope of circumferential radius p versus proper radial distance) is 
generally not far from unity except very close to a black hole horizon, where it tends to zero. 

Setting v — in (13), and using (9) and (14), one obtains the outward force needed to hold unit 
mass in a ZAM orbit: 

-Fzam = d\nV/ds r . 
This generalizes to non-equatorial ZAM orbits: its covariant form is 

(^a)zAM = u a \p vP = d a In V. (16) 

(A simple derivation is sketched in Appendix A.) 

ZAM photons emitted outward from deep in the field reach infinity redshifted by a factor V^ 1 . 
At a black hole horizon, V becomes zero, the ZAM orbits become the horizon's lightlike generators 
and the "redshifted force" VFzam becomes the surface gravity k. Equation (16) and these useful 
properties define the sense in which V = (—g tt )~i functions as the scalar potential appropriate to 
ZAM observers (see, e.g., Israel 1983). 

For definiteness, we take the angular momentum of the hole (or other central mass) to be 
positive. Then ZAM observers are dragged in the positive-^ direction (f2# > 0), at least near the 
source. The inequalities 

a > 0, (3 > 0, 7>0 (17) 

must hold, at least if the source is isolated. According to (10), these just express the conditions 
(respectively) that, as one moves out from the source, circumferences of equatorial circles get larger, 
effects of frame-dragging (as measured by fig) get progressively smaller, and that ZAM observers 
always experience inward gravitational pull. The inequalities (17) should remain valid for a larger 
class of (non-isolated) sources, in particular those for which a surrounding accretion disk or torus 
has a mass appreciably less than the central mass. 

Assuming the validity of (16), we infer at once from (9) that 

\v-\ > v+, 

showing that, at a radius r_ ph where retrograde Keplerian orbits become lightlike, prograde Keple- 
rian orbits must still exist. Thus, quite generally there is an "outer rotosphere," characterized by 
V- < —1, v+ < 1, in which retrograde Keplerian orbits (only) have become superluminal, and in 
consequence the Rayleigh criterion is reversed for counter-rotating disks and tori. 

Retrograde orbits can still extend into the outer rotosphere if they are supported, for instance 
by the pressure gradient in a thick disk. Counter- rotation in the outer parts of such a disk must, 
however, give way to co-rotation in the inner parts, since no outward radial force is available to 
hold up the equatorial inner edge. 
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Two-way structures of this kind may arise astrophysically when a retrograde inflow first impinges 
on the outlying parts of a pre-existing co-rotating accretion disk, and may survive on the order of 
an accretion time. (At least in principle, one might conceive of steady-state configurations in which 
magnetic coupling of the disk to the forward-spinning hole supplies the torque needed to drive 
the specific angular momentum I from negative to positive as the gas slowly spirals in. By itself, 
the Rayleigh criterion would forbid such a transition, since it requires \l\ to increase inwards for 
retrograde motion in the outer rotosphere. However, Archimedean buoyancy due to the denser 
inner parts of the disk can easily offset this. For compressible fluids, H0iland's criterion is the one 
that is relevant (see Blaes 1987).) 

The inner parts of disks and tori are constrained by stability and Bernoulli's law, and must stop 
well short of the corresponding (co- or counter-rotating) photon orbits. 

In thin (pressureless) disks, all orbits are Keplerian, and the innermost orbit cannot lie within 
the last stable Kepler orbit: r = r s . Stability requires that the Keplerian angular momentum \£k\ 
should increase outward. From 



one sees that \£k\ becomes infinite at both ends of the range r ph < r < oo, and hence must attain 
a minimum at a radius r s > r p h, determined from 



For the special case of the Kerr geometry, simple explicit expressions are obtainable for equatorial 
circular orbits in both pro- and retrograde cases, and are summarized in the next section and 
Appendix B. 

In the case of "thick disks" (tori), internally supported by fluid pressure, material at the inner 
equatorial edge follows a nearly Keplerian orbit subject to the condition that its binding energy, 
(1 —£k)i be positive (Kozlowski, Jaroszyriski & Abramowicz 1978; see also Appendix A). Here, the 
specific energy £ (energy per unit proper mass) is given by 



For pro- and retrograde Keplerian orbits around non-extremal holes, E+k and S-k drop from 
infinity as one moves outward from the corresponding photon orbits at r +ph and r_ ph , attain minima 
at radii r +s and r_ s respectively and finally tend to unity from below as r — > oo. It follows that 
there exist radii r e & (e = ±) such that r ep h < r e & < r ts and £ e K{r e b) = 1- Inner edges of pro- and 
retrograde thick disks cannot fall within the corresponding radii of zero binding energy r+& and r_&. 
This condition is less restrictive than that for thin disks. 

4 Circular orbits in Kerr geometry 

The formulae and results derived in the previous section for orbits on general stationary and ax- 
isymmetric fields take a rather simple explicit form in the case of the Kerr geometry. The explicit 
formulae should be a satisfactory approximation when the mass of the accreting disk or torus is 





d(f K )/dr\ r=rs =0. 



(19) 



£ = -uatfa = \ 9L/di = (1 - v 2 rHv + ptt B v). 



(20) 
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small. We collect the most interesting results here; more detail can be found in Appendix B (also 
Bardeen 1973, Lynden-Bell 1978, de Felice k Usseglio-Tomasset 1991, de Felice 1994). 

We work in standard (Boyer-Lindquist) co-ordinates for a Kerr geometry of mass m and angular 
momentum ma. It is convenient to introduce dimensionless quantities £, a* defined by 



r = m£ , a = ma*. (21) 

For equatorial circular orbits, the pro- and retrograde Keplerian 3-velocities v + and t>_ relative 
to a ZAM observer are given by 



v f — e 



m 



2 r 2 + a 2 — 2ea(mr) 2 



K r 2 -2mr + a 2 J r 2+ e am* 
with e = ±1. In terms of the dimensionless quantities (21), 

e + a 2 - 2ea^ 



V f — € 



The pro- and retrograde photon orbits v e — e have radii r ep h = Tn^ 2 ph , where 

£e P 7i(3 - C ph ) = 2ea*. 

(Another (prograde) solution is C, +P h = 0, corresponding to the singular equatorial ring r = 0, 6 = 
n / 2 inside the horizon. Although of no astrophysical relevance, this has interesting implications for 
the source structure of the analytically extended Kerr geometry (Israel 1977).) For the Schwarzschild 
case (a* = 0), Q ph = 3, and we recover the familiar result r p h = 3m for both pro- and retrograde 
photon orbits. For maximally rotating Kerr (a* = 1) we obtain ^ ph = 2, ^ +ph = 1: the retrograde 
photon circles out at radius 4m, the prograde one on the horizon itself, at radius m. 

The innermost stable Kepler orbits (pro- and retrograde) are at radii r ts = m£ 2 s , where 

U4 - (3£ e 2 s - 2)^] = 3ea*. 

This gives £± s = 6, i.e., r s = 6m for Schwarzschild. For maximal Kerr, £_ s = 3, £ +s = 1. Counter- 
rotating thin disks have inner edges beyond r_ s = 9m, far outside the outer rotosphere, which 
terminates at r_ p h = 4m. On the other hand, co-rotating thin disks can extend within the ergo- 
sphere, almost to the horizon. 

Keplerian orbits have energy per unit proper mass 

£lK = + r . (22) 

{i« 3 -3? + 2«j,)i 

Orbits of zero binding energy, 1 — Ske = 0, which arc believed to fence off the inner edges of thick 
disks (see Kozlowski et al. 1978 and Appendix A) are at radii r tb = m£ 2 b , where 

&b(2 - &{,) = ea*. 
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For Schwarzschild, = 4m. For retrograde orbits in extremal Kerr, r_6 = (3 + 2y/2)m. Prograde 
orbits never reach zero binding in extremal Kerr; evaluation of (22) gives 

1 _ S +K = 1 - = 0.423 for a, = 1, r = m, 

a well-known result (Bardeen 1970b), which gives the maximum energy extractible from accretion 
onto a Kerr black hole. 

At the radius T- V h of the retrograde photon orbit, the specific energy of the prograde Kepler 
orbit is 

£Mr-*) = ^= Jrf^X' - 3) = 2a, 

This gives a positive binding energy for (a/m) > 8^— 11 = 0.3137. If the hole is spinning faster 
than this (very moderate) value, a thick disk whose outer parts are retrograde can extend into the 
outer rotosphere provided its inner edge is co-rotating. That such hybrid objects can play more 
than a transitory role in the fuelling of quasar activity by retrograde accretion is a priori perhaps 
unlikely, but really deserves further study. 

5 Concluding remarks 

We hope that the simple treatment we have presented will contribute to the understanding of 
the general-relativistic dynamics of the inner parts of accretion disks near spinning black holes, in 
particular, the effects of frame-dragging, which, as we have seen, is very different for co-rotating 
and counter-rotating disks. Recent observational evidence suggests that retrograde accretion may 
occur sometimes or even fairly often in galactic nuclei. 

We have touched in passing on an issue which has recently had considerable exposure: Are the 
paradoxical effects inside circular photon orbits best "explained" as a reversal of centrifugal force, or 
is it preferable to suppose that kinetic energy has weight? Heuristic explanations are a subjective 
matter in which each individual is free to prescribe for himself. Simplicity (a purely subjective 
criterion) will be the deciding factor in each case. An explanation so subtle that it needs explaining 
loses its raison d'etre. We cede the last word to Richard Feynman: 

'What do I mean by understanding? Nothing deep or accurate — just to be able to see some of 
the qualitative consequences of the equations without solving them in detail." (Feynman 1947) 
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A Circular orbits in stationary axisymmetric geometry: 
covariant approach 

For any steady motion in a circle (not necessarily equatorial) centred on the symmetry axis, the 
normalized 4-velocity is 

u° = U- 1 U a , U a = ^ t) + ^ v) , (Al) 

with 

-U 2 = U a U a = g w n 2 + 2Q 9ipt + g tt . (A2) 
Axial symmetry and stationarity imply 

vPdpU = u^dpVl = 0. (A3) 

The force on a unit mass in this orbit is 

F a = u al pu = U- 2 U alf3 UP. 

Substituting from (Al) for U a \p and using (A3) and Killing's equations £ a \p = —C,/3\ a for both Killing 
vectors, we easily reduce this to 

F a = -u- 2 (u P \ a -^d a n)u p , 

or, by (A2), 

F a = U- 1 (d a U + £d a Q), (A4) 

where £ = u a ^^ is the specific angular momentum. 

A ZAM orbit is defined by £ = 0. Then U Z am = V and (A4) reduces to (16) in the text. 
The energy per unit rest-mass is 

£ = -u a z (t)a = -u a (u Ua - n^ )a ) = u + ne, (as) 

where we have used (Al). For a ZAM orbit, 

^ZAM = _ «ZAM^)a — UzAM = V. (A6) 

From (Al) and its analogue for U% AM we find 

u£ = u a ^ )a = (u a - u% AM )^ )a = (n- n B ) 9w . (A7) 

Similarly, subtracting its ZAMO analogue from (A2) and recalling Q B = —g v t/g w , 

v 2 -u 2 = 9w (n - n B ) 2 . 



n 



The relative 3-velocity v of the particle and a ZAM observer in the same spatial orbit is given 

by 

(1 - v 2 )^ = -u a u% AM = -U-^ {t) A M = V/U, (A8) 
where we have used (Al) and (A6). From (A8) and (A7), 

v = V- 1 (g vv> )kn-n B ). (A9) 
(A5) allows us to write (A4) in the alternative form 

F a = u-\d a £-nd a e). (aio) 

All these formulae of course extend at once to a continuous medium in steady rotation, with 
energy per unit rest-mass 

s = - Ua ^ t) = u + ne = v(i-v 2 )~* + n B e 

and angular momentum per unit rest-mass 

These formulae agree with (18) and (20) in the case of equatorial orbits. 

If the medium is a fluid, with local energy density /i and pressure P, the force is provided by 
the pressure gradient: 

F a = d a P/(/2 + P). 

In a "barytropic" fluid, fi depends only on P. Then F a is a gradient, and we have the usual stringent 
consequences exemplified by von Zeipel's theorem (Eddington 1926). In particular, the free surface 
(P = 0) of a barytropic fluid torus is also a surface of constant specific energy: £ = £ . If the 
surface is closed, fluid at the surface must have positive binding energy (£ < 1), since the fluid 
orbit at the outer equatorial edge is Keplerian or sub-Keplerian. (On the verge of accretion, the 
inner equatorial edge is expected to be cusped and slightly super-Keplerian, see Kozlowski et al. 
1978.) 

To conclude, let us briefly link up the present approach with the force formula (15) for equatorial 
circular orbits. 

Axi-stationarity and (A4) imply generally that F v — F t — 0. For an equatorial orbit, the 
poloidal component is zero also, and the force is purely radial. For a specific orbit we can take fl to 
be a given (position-independent) number. (The apparent dependence on V in (A4) is actually 
illusory, because it is cancelled by a part of VU.) Then (A4) reduces to the scalar relation 

17 - 1 d <U 2 ) = ^-{U 2 ) 



2U 2 ds r x 2V 2 (1 -v 2 ) ds r 

where ds r is the element of proper radial distance, and we have used (A8). 

According to (A2), U 2 depends quadratically on Q — Qb, which in turn is simply proportional 
to v according to (A9); all coefficients are purely geometrical, i.e., velocity-independent. Putting 
all this together, we arrive at 



2V 2 (l-v 2 ) ds r g vv 

1 dg vv , (i 

2g vv> ds r (1-v 2 ) 



L_ d 9w (v-v+)(v-v-) (All) 
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in agreement with (15). (All) itself identifies the quadratic roots v + , i>_ as Keplerian orbital 
velocities, for which F — 0. 



B More on equatorial Kerr orbits 

In Boyer-Lindquist co-ordinates, the Kerr equation 6 = n/2 has a metric of the form (8), i.e., 

ds 2 = rfdr 2 + p 2 {d V - Q B dt) 2 - V 2 dt 2 , 

with 

p 2 = q(r)/r, p 2 V 2 = A(r), Q B = 2ma/q{r) 

and 

q(r) = r 3 + a 2 r + 2ma 2 , A(r) = r 2 - 2mr + a 2 . 
The quantities defined in (10) take the explicit form 

a = (r 3 — ma 2 )/qr, (3 = ma(3r 2 + a 2 )/grA 5 
7 = m[(r 2 + a 2 ) 2 — 4ma 2 r]/qrA 

To compute 5, it helps to note the factorization 

[(r 2 + a 2 ) 2 — 4ma 2 r](r 3 — ma 2 ) + ma 2 (3r 2 + a 2 ) 2 = rq 2 . 

Then it easily follows that 

5 = (m/rA)2. 

Substituting these expressions into (9) produces the formulae given in Sec. 4 for the Keplerian 
velocities v±. The Keplerian energies 8±k are then found from (20). It helps to note the factorization 

{l _ v 2 ) l (e-3Z + 2ea^ 



m(£ 3 + ea*) V ^ 4 - 2^ 2 + a\ J 

in the notation of Sec. 4. 

The Keplerian angular momenta, prograde and retrograde, are 

„ / n 2x-i ^ 4 + a 2 -2ea4 

t € = pv e {l — v e ) 2 = me 



^(e 3 -3e + 2ea,)3 

The condition for marginal Keplerian stability, d(£ 2 )/dr = 0, factorizes as 

(£ 3 + ea,)(£ 4 - 6^ 2 + 8ea*£ - 3a 2 ) = 0. 

The vanishing of the second factor is most easily treated as a quadratic equation for a* in terms of 
£, and leads to the results stated in Sec. 4. 
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